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ABSTRACT 


The  phenomenon  of  cavity  resonance  Is  explained  as  an 
Interaction  between  the  flow  of  the  surrounding  liquid  and  the 
compression  and  expansion  of  the  contained  gas.  In  this  paper, 
the  general  solution  to  the  linearized  problem  of  a  symmetrical 
forebody  performing  symmetrical,  harmonic  motion,  in  an  infinite 
medium,  with  a  cavity  pressure  which  is  varying  harmonically.  Is 
presented.  The  phenomenon  of  cavity  resonance  behind  a  sta¬ 
tionary  slender  wedge  In  an  infinite  medium  is  studied  in  de¬ 
tail;  the  results  show  that  self-excited  pulsations  at  discreet 
frequencies  are  possible.  It  is  shown  that  the  resonance  fre¬ 
quencies  are  Independent  of  the  actual  size  of  the  slender 
wedge,  and  the  occurrence  of  resonance  depends  on  the  property 
of  the  contained  gas  through  the  equivalent  Mach  number  M*, 
which  is  defined  as  {pU®/p_a  where  p,  p  ,  U,  a  are  the 

density  of  the  liquid,  density  of  the  gas,  uniform  speed  of 
the  flowing  liquid,  and  the  speed  of  sound  of  the  gas  respec¬ 
tively.  The  first  four  resonance  frequencies  calculated  by 
the  present  theory  are  somewhat  higher  than  those  observed  ex¬ 
perimentally  by  Silberman  and  Song  (Reference  l)  for  a  normal 
plate  in  a  free  jet  tunnel,  although  the  trend  is  quite  similar. 
The  values  of  M*/M*,  where  the  subscript  n  denotes  the  nth  stage 

resonance,  compare  fairly  well.  The  theoretical  and  experi¬ 
mental  values  of  M*  cannot  be  compared  directly  because  of  the 
difference  in  the  forebody  shape,  they  are,  nevertheless,  of 
comparable  order  of  magnitude. 


HYDRONAUTICS,  Incorporated 

-11- 

SYMBOLS 

A  amplitude  of  harmonic  disturbance,  see  Equation  [4] 

a  speed  of  sound  of  the  liquid  In  the  uniform  stream 

00 

a  speed  of  sound  of  the  gas  contained  In  the  cavity 

C  =  +  1  Cg  ,  a  complex  constant 

P  (k/2),  P  (k:/2)  defined  by  Equations  [49]  and  [50] 

G(k/2)  defined  by  Equation  [51] 

/  2  ^ 

H'  '  Hankel  function  of  the  second  kind  and  zeroth  order 
o 

I  Bessel  function  of  nth  order,  n  =  0,  1,  2  ... 

n 

k  reduced  frequency  of  pulsation  =  oji/U 

I  cavity  length 

M  Mach  numoer  of  the  uniform  stream  =  U/a 

1 

M*  equivalent  Mach  number  -  (pU^/P^a^*) 

m  strength  of  source  distribution 

p  local  static  pressure 

p^  cavity  pressure 

p  ambient  pressure 

00 

"q  velocity  of  the  fluid  at  any  point  in  the  flow  field 

r  radial  distance  from  the  source 


t 


time 
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U  uniform  free  stream  velocity  parallel  to  the  x-axls 

u  x-component  of  the  disturbance  velocity 

V  cavity  volume 
c 

V  y-component  of  the  disturbance  velocity 

x,y  space  coordinates 

4 '  dummy  variable 

y^^  semi -thickness  of  the  body 

y^  seml-thlckness  of  the  cavity 

A  wave  length  = 

u  mass  of  gas  contained  In  the  cavity 

p  density  of  the  free  stream 

p  density  of  gas  contained  In  the  cavity 

S 

0  cavitation  number  =  (p^-p^)/5pU^ 

b  velocity  potential 

o)  frequency  of  pulsation 

The  subscript  "o"  and  "i.",in  general,  denote  the  steady 

and  unsteady  parts  of  the  designated  quantities  respec¬ 
tively. 
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(see  Equations  [49]  and  [50]) 

Cavity  Resonance  Condition  for  a  15°  Wedge  at  5  Feet 
Depth 
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ON  THE  PULSATION  OF  FINITE,  VENTILATED  CAVITIES 

I .  INTRODUCTION 

Recently,  Sllberman  and  Song  (Reference  1 )  observed  self- 
exclted  pulsations  of  ventilated  cavities  behind  normal  plates, 
symmetrical  wedges ^  and  hydrofoils  In  a  two-dimensional  free  jet. 
Later,  Song  (Reference  2)  analyzed  the  problem  by  neglecting 
the  effect  of  the  flowing  stream  and  at  the  same  time  approxim¬ 
ated  the  actual  geometry  by  a  cylindrical  cavity  enclosed  by 
an  annulus  of  quiescent  water.  The  resonance  frequency  thus  ob- 

X 

talned  Is  proportional  to  [in  (R  /r  )]^  where  R  and  r  are  the 
^  ^  o'  0  o  o 

radii  of  the  outer  and  Inner  cylinder  bounding  the  annulus  re¬ 
spectively.  The  relationship  between  the  actual  geometry  In- 
vol  ea  and  the  annulus  considered,  though  both  are  doubly  con¬ 
nected  regions,  is  not  evident;  certain  emperical  constants  have 
to  be  determined  from  the  experimental  data.  Flirthermore,  the 
result  indicates  that  as  the  outer  radius  recedes  to  infinity, 
cavities  cease  to  pulsate.  In  this  paper,  we  shall  show  that 
when  the  interaction  of  the  flowing  stream  and  the  expansion 
and  the  compression  of  the  gas  contained  in  the  cavity  is  taken 
into  account,  self-excited  pulsations  of  finite,  ventilated 
cavities  at  discreet  frequencies  is  possible  in  an  infinite 
medium. 

Starting:  from  the  linearised  potential  equation  for  har- 

* 

monic  disturbances  in  an  unsteady,  compressible  flow  together 
with  its  fundamental  solution  in  terms  of  Hankel  functions,  it 

♦We  are  ultimately  Interested  in  the  noise  generated  by  such 
pulsations;  thus,  the  assumption  of  compressible  fluid  is  mere 
appropriate  here. 
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is  first  shown  that  when  the  Mach  number  Is  small  and  the  body- 
cavity  dimension  small  compared  to  the  wave  length  of  the  dis¬ 
turbance  the  governing  equation  becomes  the  Laplace's  equation 
and  the  fundamental  solution  becomes  a  source  of  time-dependent 
strength.  For  the  problem  of  symmetrical  motions  of  the  fore¬ 
body  with  time-varying  cavity  pressuresj  a  distribution  of 
sources  with  time-varying  strength  has  the  proper  symmetry  re¬ 
quired.  By  applying  the  proper  boundary  conditions,  an  Integral 
equation  results.  To  obtain  the  solution,  further  linearization 
Is  applied;  that  Is,  the  time  dependent  part  of  the  forebody 
motion  and  of  the  cavity  pressure  are  assumed  to  be  small  com¬ 
pared  with  their  respective  mean  values.  The  solution  to  the 
steady  part  Is  given  by  Tulin  (Reference  3).  The  solution  for 
the  unsteady  part  Is  a  proper  combination  of  a  particular  solu¬ 
tion  and  a  homogeneous  solution  which  is  in  the  form  of  a 
travelling  wave,  such  that  the  Juncture  condition  Is  satisfied. 

For  the  problem  of  self-excited  pulsations  of  cavities, 
we  have  treated  a  cavitating  flow  past  a  stationary  wedge  In 
detail.  By  assuming  a  sinusoidal  cavity  pressure  variation  the 
source  distribution  can  be  obtained  from  the  method  stated  above. 
The  time-dependent  volume  of  the  cavity  is  obtained  by  a  double 
Integration  of  the  source  strength.  By  analyzing  the  gas  con¬ 
tained  In  the  cavity,  the  volume  variations  due  to  the  assumed 
cavity  pressure  variation  can  be  calculated  by  gas  laws.  These 
two  cavity  volumes,  one  calculated  by  considering  the  interaction 
of  the  streaming  flow  and  the  cavity,  the  other  calculated  by 
the  gas  lav;s,  must  be  compatible  with  each  other  in  order  to 
have  self-sustained  oscillation  of  the  cavity.  From  the  com¬ 
patibility  conditions,  the  reduced  frequencies  k  where 

Is  the  length  of  the  steady  cavity  and  U  is 


uniform  stream 
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speed)  at  which  resonance  may  occur,  and  the  equivalent  Mach 

1. 

numbers  M*  (=  (pU^/p  a  where  p  Is  the  density  of  the  uniform 

S  S 


stream,  and  p  and  a  are  the  density  and  the  speed  of  sound 
S  S 

of  the  gas  Inside  the  cavity)  corresponding  to  these  frequencies 
can  be  determined.  It  Is  found  that  the  resonance  frequencies 
are  Independent  of  the  wedge  size.  The  first  four  resonance 
frequencies  calculated  by  the  present  theory  are  higher  than 
those  observed  experimentally  by  Silberman  and  Song  (Reference  l) 
for  a  normal  plate  in  a  free  jet  tunnel,  although  the  trends  are 
quite  similar.  The  value  of  where  the  subscript  n  denotes 

the  nth  stage  resonance,  compare  fairly  well.  The  theoretical 
and  experimental  values  of  M*  cannot  be  compared  directly  be¬ 
cause  of  the  difference  in  the  forebody  shape;  they  are,  never¬ 
theless,  of  comparable  order  of  magnitude. 


HYDRONAUTICS,  Incorporated 

-4- 

II.  METHOD  OF  ANALYSIS  -  LINEARIZED  THEORY 

A  first  approximation  for  the  two-dimensional  unsteady 
cavity  flow  around  an  obstacle  In  a  uniform  Infinite  stream 
can  easily  be  obtained  by  assuming  that  the  changes  in  magni¬ 
tude  and  direction  of  the  velocity  U  of  the  undisturbed  flowj 

00 

due  to  the  presence  of  the  body,  are  small;  more  exactly,  by 
assuming  that  the  squares  and  the  higher  power  of  the  perturba¬ 
tion  velocity  can  be  neglected  when  compared  with  the  square 
of  the  uniform  stream  velocity.  The  present  analysis  is  essen¬ 
tially  an  extension  of  M.  P.  Tulin's  work  on  steady  two-dimen¬ 
sional  ca;lty  flows  (Reference  3). 

1.  Equations  ol‘  Motion 

Consider,  in  general,  the  uniform  two-dimensional 
flow  of  an  Inviscid,  compressible  fluid  of  Velocity  U  past  a 
symmetric  body  of  unit  chord  with  blunt  base  so  that  a  cavity 
is  sustained  downstream  of  the  body.  Let  the  x-axis  be  in  the 
direction  of  the  uniform  stream  and  its  origin  be  at  the  center 
of  the  base  of  the  forebody,  see  Figure  1,  Let  y^(x,t)  and  y^(x,t) 

denote  the  body  surface  and  ca /ity  boundary  respectively.  The 
pressure  In  the  cavity  p^(x,t)  is  assumed  to  be  a  given  func¬ 
tion  of  time.  Both  the  length  of  the  cavity,  i(t)  and  the 
volumeper  unit  width  of  the  cavity,  are  perm.itted  to  vary  with 
time,  the  equation  of  motion  is 

^  )  q  =  -  ^  V  p  [  1  ] 

in  which  "q  --  U  +  Vti,  <5>  Is  the  perturbation  potential,  p  and  p 
denote  the  pressure  and  density  respectively  in  the  flow  field. 
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Por  Irrotatlonal  motion  of  barotroplc  fluid,  [1]  can  be  written 
as 


d<t>  dp 

^  +  — +  p  — 


[2] 


The  perturbation  potential  <))  satisfies  the  following  equation 


M  ^ 


0^0 

dt® 


CO  O  0 

a  cjxdt 


1  d^0 

a  ®  dt^ 


0 


[3a] 


Where  M  (=  —  )  is  the  Mach  number, 
'  a  ' 

CO 


a  is  the  speed  of  sound  in 

CO 


the  undisturbed  fluid.  If  harmonic  disturbances  are  assumed, 
the  funuamental  source  pulse  solution  to  [3]  Is  given  by  Bis- 
pllnghoff,  Ashley,  and  Halfman  (Reference 


^Vl-M  ^ 


exp  1 


M 


a)t+ 


CD 


1-M 


(x-g.) 


H 


(O 


— - V  (x-g)^"+  (3-M/)ry-Ti)^ 

a(l-M  V  . 

uO 


in  which  (g,q)  is  the  source  point,  and  A  is  the  amplitude  of 
the  harmonic  disturbance.  For  M  «  1,  which  Is  of  interest 

CO 

here,  the;  above  expression  becomes 


0  (x,y,t) 


lAiliHl 

4 


o 


in  which  r 


+  (y~q)^-  it  is  known  that  (see  for 


example  Jahnke  and  Eiide,  Reference  5) 
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'  If 


< 

„i(r/A-7r/4) 

irr 


for  ^  «  1 


[4a] 


[4b] 


If  the  body-cavity  dimension  is  much  smaller  than  the  wave 
a 

length  "K  =  ,  In  the  near  field  of  the  body,  the  fundamental 

solution  assumes  the  form  of  a  source  with  time  varying  strength: 


<t>  =  ^ ^  ^  e^^^  in  (r/A)  [5] 

which  satisfies  the  governing  equation  for  an  incompressible 
fluid: 

<t>  -  0  [6] 

(Equation  [6]  is  obtainable  from  [3a]  by  the  assumption  of 

0  (r^/A^)  =  0  (M^r/A)  which  are  negligible  compared  to  1.] 

As  r  increases,  ultimately  the  form  of  <t>  in  [5]  must  be  replaced 
by  the  expression  [4b]  which  vanishes  as  r  This  behavior  in 
the  far  field  is  due  to  the  effect  of  compressibility  which  be¬ 
comes  dominant  on  a  scale  larger  than  the  wave  length  of  the 
pressure  disturbance.  Therefore,  tlie  boundary  con.  ition  that 
the  pressure  obtained  by  an  incompressible  analysis  must  be  zero 
at  infinity  as  used  by  some  investigators  is  not  only  unnecessary 
but  in  fact  erroneous. 

The  pressure  integral  in  [?]  may  also  be  replaced  by  p^ p 
and  the  pressure  at  any  instant  may  be  obtained  by  integrating 
[2] 
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Pc°~P  ^  2  I  A. 

1  ..2  U  Sx  U  St 
2p  U  ' 


[7] 


If  the  same  restrictions  are  placed  on  the  gas  inside 
the  cavity,  l.e.  ca£/a  «  1  M  (=  U/a  )  «  1,  where  a  is  the 

s  s  s  s 

speed  of  sound  of  the  gas.  We  have,  essentially,  an  incompres¬ 
sible  unsteady  flow  problem  in  the  near  field  of  the  body-cavity 
region . 

2.  Boundary  and  Closure  Conditions 

The  condition  at  the  boundary  of  a  typical  body 
states  simply  that,  over  its  surface,  the  normal  component  of 

the  fluid  velocity  ^  is  fixed  by  the  body's  motion.  If  the 

equation  of  the  surface  of  a  body  moving  in  a  time-dependent 
fashion  is 


(x,t)  -  0  -  1  <  X  <  0  [6] 

then  the  linearized  boundary  condition  to  be  satisfied  along 
the  x-axls  is 


v_ 

U 


i  A. 
u  at 


y^,  (x>t) 


l<x<0,  y  =  0 


[9] 


Where  v  is  the  velocity  component  in  y-direction.  On  the  cavity, 
the  pressure,  p^(t),  is  specified,  or  equivalently,  the  cavita¬ 
tion  number  a  is  specified.  From  [7],  it  follows  that 


P  -P„(t) 

00  C 


-  aft) 


1  A. 

u  at 


ip 


2 

u 


c  <  X  <  i,  y  =  0 


[10] 
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The  condition  that  the  cavity  be  closed  at  every  Instant 
can  be  derived  from  the  usual  klnematlcal  requirement  that  at 
a  bounding  free  surface  the  motion  of  the  surface  must  coincide 
at  any  time  with  the  motion  of  those  fluid  particles  that  happen 
to  be  at  the  surface  at  this  time;  that  Is, 


hi 

U 


d4  =  0 


[11] 


We  further  specify  that  the  motion  of  the  forebody  be  symmetric, 
with  respect  to  the  x-axis. 


3.  General  Solutions 


The  linearized  mathematical  problem  may  be  stated  as 
follows:  To  find  a  harmonic  function  <t>(x,yjt),  symmetric  with 
respect  to  the  x-axis,  whose  gradient  in  the  limit  vanishes 
everywhere  on  a  circle  of  sufficiently  large  radius  about  the 
origin  which  satisfies  the  mixed  boundary  conditions  [9J  and 
[10]  and  the  closure  condition  [11].  A  distribution  of  unsteady 
sources  of  strength  K,'..  t)  along  the-  x-axis  I’or  -  1  <  x  <  -^  pro¬ 
duces  a  harmonic  function  with  the  proper  symmetry.  Therefore, 
we  have 


♦  (x.y.t)  ^  — 


m  (^,t)  £n 


-1 


with  velocity  component: 


1?] 


u  (x,y;t) 


b<t>  1 

dx  27r 


t) 


(x-?)  ae  l13J 
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V  {x,y;t)  =  ^ 


1 

2Tr 


-9- 

r 

-1 


-  _ 

(x-O^+y^ 


y 


[14] 


Here  and  subsequently,  £  =  It  can  be  shown  that  at  the 

x-axis 


u  (x,0,t)  =  — 


x-^ 


di 


[15] 


-1 


/y  (x,  t)  -  l<x<0 


V  (x,0,t)  =  •|■m(x,t)  -  U 


l^x  U  ht 


16] 


[y  (x,t)  0<x<i 


The  mixed  boundary  condition  on  the  x-axis  will  be  satisfied  if 


m(x,t)  =  2U 


+  i  A 

dx  U  dt 


yb(x,t) 


-  1  <  X  <  0  117] 


and 


^  ^  I  •^n(x-0  d|  -  TrU0(t)  C  <  X  <  £  [lb] 

-1 


Decomposing  the  source  distribution  into  steady  and  unsteady 
parts,  and  using  the  subscripts  "o"  and  "1"  to  denote  the  steady 
and  the  unsteady  part  respectively 


m(x,t)  m  (x,t)  =  m  (x)  +  m  (x,t). 


m 


b,i 


^  ^  ^  «1  1  <  X  <  0  [19] 

’  ■*  b ,  o 


m 


m(x,t)  =m  (x,t)  -  m  (x)  +  m  (x,t)  «1  0  <  x  <  i  [20] 


c,o'  '  C,1 


G,  O 
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and 


«  1  0  <  X  <  i 


[21] 


[l8]  becomes 


0 


TrUa  = 
o 


and 


ttUo  (t )  = 


m  (^) 

x-e 


d4  + 


m  (^) 


c,  o 


d^,  0  <  X  <  ,  [22] 


0 


^  +  1  A. 

dx  U  dt 


J 

-1 


0  ^ 

%^i(?.t)in  x-^  +J  m^^^(^,t)£n 

0 


x-4 


0  <  X  <  i 


[2 


in  which  Z  is  the  length  of  the  steady  cavity  and  from  [17]> 
b,o  dx 


,11.  1  1.1 

\\,x  ■  [dx  ■■  U  dt 


'  b,  1 


[24] 


[25] 


and  it  is  assumed 


b,  o 


«  1. 


The  general  solution  for  [22]  is  given  by  Tulin  (Reference  3) 


m  (x)  = 
c .  o 


Tr‘ 


0 


x-^ 


ttUo  -2U 
o 


r 


dx 


d4 


)  [26] 


-1 
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Por  the  solution  of  [23],  we  first  rewrite  the  equation  as 


de  =  7rUa^(t) 


!  "b,x 


[23a] 


It  Is  noticed  that  any  arbitrary  travelling  wave  solution  can 
be  added  to  the  particular  solution  vjlthout  disturbing  the 
function  given  on  the  right  hand  side  of  the  equation,  i.e, 
this  equation  admits  homogeneous  solutions  of  the  following 
form 


m^^  (x,t)  •-  - 

c,i '  ' 


ir^  vxu-x 


(x-0 


f’  (x-Ut)  d^ 


which  is  the  solution  of  the  following  integral  equation 


m*^  (?,t)  Zn  x-^i  d^  f  (x-Ut). 

c ,  1 


A  particular  solution  of  [23a]  is 


(x,t)  = 


'x(i-x) 


u/Ua^d)-  ^  /  '"^^,(^,t)in|x-^id^  dt 


where  B(t) 
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which  may  be  verified  by  substitution  of  [28]  and  [29]  Into 
[23a]  and  noting  that 

i 

^  d|  =  0. 


By  contour  integration,  it  can  be  shown  that 


Jin  |x-g|  dj 


Zn 


I 

4 


[30] 


we  note  then  the  particular  solution  becomes  singular  at  i  =  4. 
The  solution  given  here  then  is  valid  only  for  i  7^  4,  This 
point  will  be  discussed  further  in  a  latei>  section. 

The  unsteady  source  strength  is  obtained  by  a  ccmbina- 
tlon  of  the  homogeneous  and  the  particular  solutions,  then  re¬ 
quiring  that  the  juncture  condition  holds.  This,  however,  does 
not  lead  to  a  unique  solution  in  general  since  there  is  an 
arbitrary  function  involved.  For  simple  harmonic  oscillations 
of  the  body  and/or  of  the  cavity  pressure,  we  may  assume 


P' (x-Ut) 


Re 


;;  ia)(t-x/u) 


.31] 


where  C  is  a  complex  constant 

c  -  +  i  ,  [32] 

and  Re  denotes  the  real  part.  Then  the  juncture  condition  is 
sufficient  to  determine  the  constants  and  .  The  general 
solution  of  [23]  is 
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m 


c,i 


TT^  '\/x(i-x)  (_ 


£ 


[33J 


The  cavity  shape  may  now  be  found; 


X 

r 


jjx,t)  =  y,(0,t)  +  I  d^l 

-i  -b 


[3^] 


The  closure  condition  [11]  becomes,  for  i  <  i 


y^(0,t)  +  ^ 


f 


-0  ’  -* 


for  £  >  £ 


2u  j  "’c,!^^’^ 
0 


[  36] 


since 


h 

—  hi  "c,o(^) 


137] 


The  volume  of  the  cavity  per  unit  span  is 


V  (t) 
c '  ' 


/ 

■  ij  <*’■  / 


/ 

0 


“c.o'?)  +  '"cf?-''  ■  ¥> 


0 


38] 
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III.  CAVITY  RESONANCE 

To  analyze  the  cavity  resonance  problem,  we  focus  our 
attention  to  the  problem  of  a  stationary  wedge  with  a  trailing 
cavity  maintained  by  ventilation.  It  Is  assumed  that  the 
cavity  pressure  Is  oscillating  about  a  mean  pressure;  the  reso¬ 
nance  condition  corresponds  to  this  case  when  such  oscillations 
are  possible  with  the  mass  of  gas  within  the  cavity  being  con¬ 
stant  . 

Equation  [38],  together  with  Equations  [26]  and  [33]^ 
gives  the  response  of  the  cavity  volume  to  the  time-varying 
cavity  pressure.  At  the  resonance  condition,  this  volume  varia¬ 
tion  must  be  compatible  with  those  produced  by  the  gas  Inside 
the  cavity. 

Let  p  be  the  constant  mass  ano  p  (t)  be  the  time-dependent 

Ej 

density  of  the  gas  in  the  cavity.  The  volume  of  the  ca/ity  is 


V  (t) 
c '  ' 


..F— .. 


The  variation  of  the  volume,  V  (t),  due  to  a  small  variation 

C  ,  1 


of  pressure  is 


(t)  =  - 


[39] 


Assuming  that  the  pressure  variations 
change  of  the  density  with  respect  to 

dp, 


sidered  isentropic,  then 


S 


dp. 


are  not  too  rapid,  the 
the  pressure  may  be  con- 

Equation  [39]  becomes. 


when  the  'ariation  in  cavity  pressure  i.,  written  in  terms  of 
the  cavitation  number 
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V 


(t) 


V 

c. 


o 


1 

2 


pir 


p  3 

g 


a,(t)  I  a^(t) 


[40] 


in  which  the  equivalent  Mach  number  M*  is  defined  as 


M*^ 


^  tt2 

HP 


pa' 

^g  g 


[41] 


Equating  the  volume  V  (t)  calculated  by  [40]  to  that  obtained 

C  ,  1 

by  [3hJj  the  condition  at  which  resonance  occurs  can  be  deter¬ 
mined. 

Let  the  wedge  profile  be  described  by 

yb(x)  (1+x), 


and  the  pressure  in  the  cavity  by 


o(t)  ^  0  (l  +  '"  sin  cot)  ['t3] 

'  o 

where  co^  and  i.u  are  the  amplitude  and  frequency  of  the  fluctua¬ 
tions  of  the  ca  1  ty  pressure  re.specti -  e] y ,  Then 


c(t) 


iMl 

la 


0  cos  ait  . 
o 


[44] 


To  apply  the  juncture  condition  the  integral  in  !3i]  at  x  -  0 

*  p  0 

can  be  evaluated  by  using  the  angular  transformation  ^  ^  sin  — 
and  the  Jacob's  expansion 


^i  k/2  cos  0 


J  (k/2)  +  ?  ''  i"  J  (k/?)  cos  n  0  . 
o  '  n  '  ' 

nTi 


We  obtain  for  the  unsteady  source  strength 
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m 

Cjl 


o 

k 


r 


in  Y  ~\fx(£-x) 


4 

i 


i7r(J  ^2) 


'Vx'  (i-x'  ) 


cos  cot  + 


„  "V  lcJo(t-X'/U)  -  ar'r-n 
Re  C  e  '  ^  ■'  ax '  >  [  ^5] 


X  -  X' 


0 


where  C  =  C  +  1  C 
1  2 

(J^  sin  ^  f)  +  ^  f  ” 

k  =  ^  ,  the  reduced  I'requency, 

and  the  arcuments  of  the  Bessel  Functions  and  are  k/2. 

To  apply  the  closure  condition  to  find  the  relationship 
beoween  the  unsteai.iy  ca /ity  xoneth  i(t)  =  i^+  i^(t)  and  the  un- 

steauy  cavity  pressure,  we  first  use  37]  in  351  to  obtain 


!|. 


c ,  o 


(O  ae 


.  Azi 


) 


Con. istent  with  the  assumptions  made  alrea 

- 1 -  «  1.  The  first  inte^vral  can  be 

o 

all  x's  in  the  m  expression  ^26]  by  i 
c ,  o  o 

_  1 

part  (i^-x)  ^  .  The  result  is 


dy,  we  assume 
eva-vuated  by  replacing 
except  the  singular 


const.  V^ij^(t) 
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The  second  Integral  is  clearly  of  order  g.  Equation  [47]  gives 
the  result  that 


|^i(t)|~  0  (g^)  , 


and  it  may  be  neglected  in  our  present  theory. 

In  view  of  the  above  result,  to  obtain  the  additional 
cavity  volume  due  to  the  unsteady  pressure,  we  may  replace  I  by 
in  the  expression  [45]  for  ^  including  that  in  the  reduced 

frequency  k  and  in  the  upper  limit  of  the  integral  for  the  cavity 
volume : 


V 

c,i 


(t) 


u 


)  cie  . 


This  Integration  can  be  carried  out  once  the  angular  transform:- 
tion  ^  -  sin®  —  and  the  Jacobi's  expansions  are  used.  The  re¬ 
sult  is 


V 

c,i 


(t) 


GO  Umi* 

o 


k‘ 


(k/2) 


sin  a>t 


(k/2  )  cos  cat  " 


[4bJ 


in  wluch 


F^(k/2)  =  cos  ^  -  1  +  G  (k/2)  [ sin  ^  cos  ^  J  [49] 

F^(k/2)  =  sin  I  -  G  (|)  [J^  cos  |  sin  |  ]  [50] 

G(k/2) .  [J  ^  y "  ■'7 ' 


and  all  the  arguments  of  the  Bessel  functions  arc  k/2.  The 
additional  volume  given  by  [ 4o]  becomes 
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V 


V 


e  a  sin  cut 
2  o 

0,0 


[52] 


For  small  values  of  a^,  Tulin's  result  (Reference  37)  gives 


V 


Cj  0 


TrJi  ^ 

-T-  “o 


[53] 


Combining  [53]  with  [4b],  and  comparing  with  [52],  we  require 
that 


k 

F  (tr)  =  ^ 
2  2  ^ 


*2 


a  in  ^ 
o  4 


-  —  F,  if) 

io  ^ 


[54] 


l  55] 


The  subscript  n  denotes  the  nth  1^00^  oi'  F  (k/2)  for  which 


^  (V^) 


(-jp) 


>  0  . 


This  means  that  of  >  4,  resonance  may  occur  when  (k/2)  >  C. 

Whereas  for  ,  the  resonance  conuition  becomes  Fj^(k/2)  <  C. 

For  cavity  lengths  in  the  neighborhood  of  four  chord  ien^.ths  of 
the  wedge,  M*  becomes  arbitrarily  large.  Since  the  value  of 

—  F,  (■^)  is  highly  damped  as  k  becomes  lar.  e,  only  high  modes 
k^  ^ 

of  resonance  may  occur  in  the  vicinity  of  4 .  This  result 

is  shown  in  Figure  3>  vjhlch  ivill  be  discussed  presently. 

From  [50]  and  [5l]>  it  is  seen  that 
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(k/2)-^k^ 

\  as  k  0 

—  P^  (k/2)  -  3 
k^ 

which  means  that  for  the  case  >  4,  k  =  0  is  an  admissible 
root.  But  this  corresponds  to  cd  =  0  and  V  =0,  v;hlch  re- 

C  ^  1 

duces  to  the  steady  case.  It  is  seen  from  [9^]  that  the  reso¬ 
nance  frequencies  are  independent  of  the  VMedge  size.  The  func¬ 
tion  ^  P^  (l</2)  together  v\/lth  —  P^  (^)  are  plotted  in  Pigure 

k^ 

2.  It  is  seen  that  for  Ji  >  A  the  first  non-zero  root  as  v.'ell 

o 

as  the  third,  the  fifth,  and  so  on  of  P^  are  inaumissible  be¬ 
cause  the  corresponding  value  of  P^  is  negative.  The  inoumissi ble 
roots  for  the  case  of  £  >  ''  become  the  admissible  roots  for 

the  case  of  £  <  .  The  first  four  admissible  roots  of  P  (k/2) 
for  both  £  >  ■(  and  ^  <  l  anci  the  correspond! n.';  /aj.ues  of 

CO 

/k\  ,  ,  mt-  ^ 

—  F  (— )  are  [liven  in  TabLc  I, 
k^  "  ^ 

Silberman  and  Song  ''Reference  l)  have  presentee  experi¬ 
mental  data  of  pulsating  cavities  behind  n  l/O”  normal  plate 
in  a  free  Jet  tunnel.  They  have  observed  cavity  pulsition  which 
are  characterized  by  the  number  of  -waves  appeariig  on  the  ca.ity. 
They  have  designated  the  pulsation  with  n  waves  a.,  the  nth  sta^e 
pulsation.  It  is  conjectured  iiere  that  tne  present  nth  k  would 
give  a  similar  cavity  shape.  Since  the  1  n'lefin.ite  integral  for 
the  cavity  shape  is  extremely  difficult  to  evaLuate,  It  has  not 
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been  attempted  here.  Prom  these  data,  the  reduced  frequency  and 
the  equivalent  Mach  number  corresponding:,  to  each  of  the  stages 
may  be  calculated.  In  calculating  the  density  of  the  contained 
air  at  the  cavity  pressure,  we  have  assumed  that  the  Isothermal 
equation  of  state  holds.  These  values  are  also  given  In  Table  I. 
The  reduced  resonance  frequencies  are  given  as  a  range  because 
the  experimentally  observed  cavity  length  varied.  The  values 
of  have  also  been  calculated  and  presented  In  the  Table  I. 

Since  In  the  theory  the  cavity  length  Is  measured  with 
respect  to  the  chord  of  the  wedge,  the  experimental  results  ob¬ 
tained  from  a  1/b"  normal  plate.  In  a  strict  sense,  cannot  be 
compared  with  the  theoretical  results.  However,  it  Is  of  in¬ 
terest  to  see  the  orders  of  magnitude  of  the  values  obtained. 

With  this  In  mind,  we  make  the  follov;ing  comparison.  The 
theoretically  predicted  resonance  frequencies  are,  in  general, 
higher  than  those  observed  experimentally.  The  values  of 
MVM,*  seem  to  be  in  fairly  good  agreement.  The  theoretical 

results  are  for  an  infinite  medium.  The  effect  of  the  free 
surfaces  bounding  the  jet  in  the  case  of  the  experiment  would 
probably  alter  the  theoretical  results  to  some  extent.  Prom 
the  results  obtained,  however,  it  is  seen  that  the  presence  of 
the  free  jet  surface  is  not  a  necessary  condition  for  the  occur- 
rance  of  self-excited  cavity  pulsation. 

In  terms  of  numerical  examples,  a  l'J)°-wedge  with  a 
trailing  cavity  at  a  mean  cavitation  number  of  0.07>  the  equiva¬ 
lent  Mach  number  of  the  first  stage  pulsation  is  1.25  as  compared 
to  the  experimental  value  of  2.26  for  a  normal  plate  in  a  free 
Jet  with  a  cavitation  number  of  0.0735*  It  is  interesting  to 
note  that  these  are  of  the  same  order  of  magnitude  in  spite  of 
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the  different  forebody  confii£,uratlons  and  the  free  surface  con¬ 
ditions  present  In  the  experimental  case.  In  the  presence  of 
free  surfaces.  It  Is  known  that  the  cavity  length  Is  shorter 
than  that  In  an  Infinite  medium  (Reference  6).  It  Is  con¬ 
jectured  here  that  the  effect  of  the  free  surface  Is  to  Increase 
the  value  of  the  equivalent  Mach  number  because  of  the  decrease 
In  cavity  length.  This  15°-wedge  ;\ilth  =  0.07^  may  experience 
the  first  stage  self-induced  pulsation  at  a  forward  speed  of 
50  ft/sec. 

In  Figure  3)  we  have  shown  the  locus  of  the  points  for 
which  resonance  of  the  first  three  stages  are  possible  for  a 
15°-wedge  travelling  at  5  feet  below  the  free  surface.  (in 
this  example,  we  assume  the  presence  of  the  free  surface  does 
not  alter  the  results  obtained  above.)  Two  families  of  loci 
appear,  one  for  <  4,  one  for  .2^  >  4 .  Cavity  resonance  is  only 


possible  when  the  values  of  —  and  U  correspond  to  a  point  on 


one  of  these  curves.  It  Is  interesting  to  note  that  the  family 

of  loci  for  i  >  ^  exhibit  a  minimum  speed  below  which  no  re- 
o 

sonance  corresponding  to  that  state  or  lower  may  occur.  We  may 
trace  the  history  of  the  unsteady  motion  of  a  cavity  behino  a 
15°-wedge  as  follows.  Suppose  P^/P^  =  0.75.  As  speed  increases 

from  zero,  first  stage  oscillation  may  be  encountered  at  30  ft/ 
sec.  As  the  speed  is  increased  further  the  flow  becomes  stable 
until  at  50  ft/sec,  second  stage  oscillation  appears.  Then  as 
the  speed  is  increased  still  further,  higher  and  higher  modes 
of  oscillation  appear.  At  60  ft/sec,  which  corresponds  to  the 
critical  length  =  4,  no  oscillations  may  appear.  Between 

U  =  60  and  o6  ft/sec,  the  mode  of  oscillation  starts  high  and 
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decreases  to  n  =  1.  Beyond  65  ft/sec,  oscillations  of  the  cavity 
are  again  not  likely.  It  is  in  the  range  of  speeds  from  50  to 
66  ft/sec  that  self-excited  oscillations  would  seem  most  likely 
to  appear.  The  results  presented  in  this  Figure  shovj  that 
cavity  resonance  can  only  occur  for  a  finite  range  of  speeds. 
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IV.  SUMMARY  AND  CONCLUSIONS 

The  theory  presented  In  tW  c  report  gl'/es  a  simple  method 
for  determining  approximately  the  characteristics  of  t;\/o- 
dlmensional  unsteady  cavity  flows  about  slender  symmetrical 
bodies  in  a  uniform  infinite  stream.  The  present  paper  is  in¬ 
tended  to  throiv  some  light  on  the  physical  mechanism  of  the 
pulsation  of  the  cavities.  The  account  is  mainly  theoretical, 
but  a  brief  reference  is  made  to  some  recent  experiments  on 
the  instability  of  ventilated  calties.  The  compressibility  of 
the  gas  inside  the  cavity  is  pointed  out  to  be  an  essential 
factor  in  this  aspect  of  cavitation.  The  case  of  ventilated 
cavity  flovj  about  a  thin  vjedge  has  been  discussed  in  detail. 

The  important  results  obtained  in  this  study  may  be  sum¬ 
marized  as  follov\fs: 

1.  Pulsation  of  finite,  .-entilated  cavity  is  possible 
in  an  Infinite  medium  under  ordinary  conditions  on  the  speed 
and  cavitation  number.  Within  the  framev;ork  of  the  analysis 
made,  the  existence  of  a  free  surface  is  not  a  prerec.uisite 
for  self-exclteu  pulsation. 


2.  The  occurrance  of  resonance  depends  on  the  equivalent 
Mach  number  M*  which  Involves  not  only  the  property  of  the  flo..'- 
ing  liquid  but  also  that  of  the  contained  gas  in  the  cavity: 

i 

M*  =  (p/p^)'"  U/a  _  .  The  reduced  resonance  frequencies, 

h 


k  -  cjj£^/U  ,  are  Infinitely  many  and  they  are  discreet. 

3.  The  resonance  frequencies  calculated  for  a  ca  ity 
created  behind  a  tv;o-dimensional  ueoge  is  iiuepenuent  of  the 
wedge  slse.  The  equivalent  Mach  nunoer,  however,  depends  on  the 
length  of  the  ca  /ity. 
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4.  Calculations  made  for  stationary  wedges  show  that 
there  exists  a  critical  cavity  length  of  =  4.  The  resonance 

frequencies  and  the  e^iUlvalent  Mach  numbers  for  cases  In  which 

Z  >4  are  quite  different  from  those  for  I  <4.  Furthermore, 
o  ^  o 

the  theory  shows  that  cavity  resonance  can  only  occur  within  a 

finite  speed  range. 

5.  It  Is  conjectured  that  the  nth  resonance  frequency 
would  give  a  cavity  shape  with  n  waves  similar  to  those  observed 
by  Sllberman  and  Song  (Reference  l). 

6.  The  calculated  resonance  frequencies  and  equivalent 
Mach  numbers  are  qualitatively  comparable  to  those  observed 
experimentally  by  Sllberman  and  Song  (l)  In  a  free  jet  tunnel. 
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FIGURE  I-  SCHEMATIC  OF  UNSTEADY  CAVITATION  FLOW  PASTA  THIN  BODY 


3AT 
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TABLE  I.  THEORETICAL  AND  EXPERIMENTAL  RESULTS 
(Experiments  of  Sllberman  and  Song) 


n 

1 

2 

3 

4 

^  >  4 

k  ° 

theory 

8.4 

14.8 

21.5 

27.4 

o 

A 

5.5 

11.8 

18 

24.2 

kexp(^  "  plate) 

4. 5-6.1 

9.9-13.0 

14.2-19 

17.1-22.9 

i  >  4 

k  i  <  4 

o 

0.192 

0.097 

0.055 

0.044 

-1.11 

-0.32 

-0.16 

-0.095 

rjT  (theory) 

1  jg  <  4 

o 

1 

0.505 

0.29 

0.23 

1 

0.536 

0.36 

0.30 

1" 

M  exp(^  normal  plate) 

2.26 

1.25 

o.biS 

0.631 

M* 

(exp) 

1 

0.55 

0.36 

0.28 

i 

i 
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